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Cosmic ray (CR) propagation in the Galaxy and general- 
ly in the cosmic plasma is usually considered in the duffu- 
sion approximation. The diffusion is regarded to result 
from CR scattering due to their interaction with a magnetic 
and an electric field. In most cases the fields are assumed 
to be given. Meanwhile, in the Galaxy the CR energy densi- 
ty 'vfj. ~I e?/cm, i.e. it is comparable with the energy densi- 
ties of the magnetic field and turbulent motions in the in- 
terstellar gas. Therefore, for the Galaxy and for a number 
of other objects it becomes necessary to take into account 
the influence of CR on the gas dynamics and on the magnetic 
fields in this gas; see /I ,2/. The simplest way to do this 
is to use the hydrodynamic approximation, but this is pos- 
sible only on scales greatly exceeding the CR free path A 
before scattering and only for times larger than A/v csA/c. 
One should thus obtain corresponding MHD equations and es- 
tablish the limits of their applicability. 

What does the CR pressure affect ? The mean force with which 
OS act on the gas is usually written as ^F> = -vP cr , where 
P cr is the CR pressure. But what is the pressure in our 
case where particles do not collide? It is clear that the 
forqe with which CR act on the background plasma is equal 
to F a c" 1, ( Her ) » where is the current in plas- 

ma and Hc r = B cr is a magnetic Afield created by the CR 
current y cr . For quasistationary fields 

(I) 

Here wejtiave integratajL overman ensemble of random magnetic 
fields H = > + s Hi, ^SH^> = 0 a &d, accordingly, j* s 

<T> +|f , <£'BT er > 2 * * * 6 = calculate u r , one should 

kn8ft the CR ^distribution function f(t ,r,3). It obeys the 
equation (we assume that the electric field g = 0) . 

( 2 ) 

Assuming f = ^ f>+ < of ,<S f > = 0 and using the standard pro- 

cedure of the quasilinear approximation, we find equations 
for ^ f> and 9f. Integrating the equation for S f in the 

approximation of "nonmagnetized" CR with a charge Ze (this 

means that in calculation of the interaction between a par- 
ticle and a single magnetic field inhomogeneity, particle 
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motion may be considered rectilinear, i.e. the gyroradius 
r H »L, where L is the correlation scale of the field Si?*} 
we obtain t d * 

= (?- (3) 


In the diffusion approximation ^ f> has the form (for 
more details see /3/) . 

?)>=£- =^ ; i V(f., /V = j>,>f„, ^cr = P 3 r^fo, (4 > 
where f (t,?,p) is an isotropic part of the function *f>, N 
is the total GR concentration, the diffusion tensor is eau- 
al to (here t = <K> 

O S 

D i .= v Af3(d + AVr, 2 )J- 1 7s ij 4(A)H;L.* ^ 

The free path along the f ield^ff > is . 

g&0 

A = 2 ^f> 2 r Zi z e z L 't(y) ] _1 

We assume that 0 

(?,)> = f ^(IKj - X(i) *-£ ), (7) 

vJhere x” = xt-X 2 , the functions ^ and 'V’j are related by 
the condition v^H~ = 0. 

Using formula! (2 )-(7) , we^find 

*F>=-vP cri ? c( , = ftr £ 0 (*. ?. p), ( 8 ) 

which is the answer to our question. No t^ that for i^»rn the 
main contribution into the total force cF> in the direction 
longitudinal with respect to the regular field is made by 
the fluctuation i nt era cti on - ^ f cr x /4/, and in 

the perpendicular direction - by the interaction - ■I <£ l s, cr> x 
(more precisely, by its part connected with the Hall 
current of OR). We have, in fact, extended the problems of 
momentum exchange between a relativistic charged particle 
and a "magnetic cloud” (see /5/) to the case of GR propaga- 
tion in a turbulent medium. 

MHD equations for interstellar medium with an account of the 
action of GR. We proceed from the hydrodynamic desc-rifrtlcm 
of back-ground plasma motions and from the description of 
GR by means of the kinetic equation. The equation of back- 
ground plasma motion in the MHD approximation has the form 

^ ^ j Vy (lj * (9) 

here H = H + H cri , £ = - \ x F? is an electric field in a 

medium. We assume that in plasma there exists an excessive 
charge density -ZeN which compensates the GR charge. 

. The averaged kinetic equation for OR now has the form 


^ *(jrV)4> + te < g’^.> + 2e (f * * 2e<(f *5«|p>=0(IO) 
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In the diffusion approximation with an account of the mo- 
tion of the medium the distribution function <f> differs 
from (4) by an additional term , i> (af-g ) k Mg . Multiplying (10) 

-s» 5 Tt*- ' ^ p 

by p and integrating over d^p , we nave 

If - «P. r - 0- (II) 

Here we elimitate terms quadratic in the small parameters 
-u/tr » A /H (R is the characteristic scale of the change in 
the OR concentration) . The quantities vr cr and P cr are de- 
fined by formulae (4), (8). Note that in hydrodynamics the 
density vr of the internal gas density and the gas pressure 
p are defined in the frame of reference moving at a mean 
mass velocity iT of the medium. In the case of OR it is mo- 
re convenient to use a laboratory frame of reference. If the 
quadratic terms (u/v)^ are neglected, both the definitions 
of the quantities ms and p are equivalent. 

Summing up equation (II) and equation (9) averaged 
over the field fluctuations, we derive 

= P« r + *?>-»• (12) 

^ In tne quasistationary approximation the regular field 
•^H> is determined from the equation 


il (*£,.> + (IJ) 

Substituting (13) into (12) and neglecting the small terms 
*5 H c* H, ( + p cr ) , we have 

$ V V + -fa fa* *#>)**$>. (I4) 

We find the equation for ovT cf by multiplying (10) by 
E and integrating over d-^p 

l^r r - v £>v v; r -t (oT^)^ r + % r ) (*5) 

The same result can be obtained simpler, by using the equa- 
tion for f 0 in the diffusion approximation:|^ , -v.D t ..v.-} t ,+^ 7 )/,-^l 4 ^ 

The rest of the MHD equation, except (14) , (15) are ^ 
obvious: ^ 

+ = ?o£ V^JP> 0 . (16) 

Equations (14) -(16) form a system of MHD equations for inter- 
stellar medium with an account of the action of OR, They must 
be supplemented with the equation of state for the backgro- 
und plasma and with the relations connecting p and -w; 


(see /6/) » For ultrarelativistic OR p 


cr 






'cr 


In the approximation used here, CR are an additional 
gas component with a large internal energy, but with a neg- 
ligible mass density. The small-scale random field S H H, 
which provides OR scattering, enters only the tensor ft? The 
field <5 H induced by the fluctuations 'Sf must be small: 


cr 
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SH cr otherwise CR scattering cannot be considered 

in the test particle approximation* This condition is fulfil 
led if ( Sc,, is the degree of OR anisotropy) 

? ct .S cr H 2 Att- d7) 

is is assumed that in the system there are neither kinetic 
instabilities induced by OR (see /?/) nor a statistical 
acceleration of particles. 

Various dissipative processes maybe included into the 
equations if necessary. Note that OR scattering leads to 
viscosity and to large-scale conductivity (along the field 

K • f 5 WM. > **, - A f„. < I8 > 

For examples of the solution of the above MHD equations 
see /2,8— 10/). 
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